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/3-decay properties of even and odd-A neutron-rich Ge, Se, Kr, Sr, Ru, and Pd isotopes involved 
in the astrophysical rapid neutron capture process are studied within a deformed proton-neutron 
quasiparticle random-phase approximation. The underlying mean field is described self-consistently 
from deformed Skyrme Hartree-Fock calculations with pairing correlations. Residual interactions 
in the particle-hole and particle-particle channels are also included in the formalism. The isotopic 
evolution of the various nuclear equilibrium shapes and the corresponding charge radii are investi¬ 
gated in all the isotopic chains. The energy distributions of the Gamow-Teller strength as well as 
the /3-decay half-lives are discussed and compared with the available experimental information. It 
is shown that nuclear deformation plays a significant role in the description of the decay properties 
in this mass region. Reliable predictions of the strength distributions are essential to evaluate decay 
rates in astrophysical scenarios. 
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I. INTRODUCTION 

The rapid structural changes occurring in the ground 
state and low-lying collective excited states of neutron- 
rich nuclei in the mass region A = 80 — 128 have been 
extensively studied both theoretical and experimentally 
(see e.g. [l|, 0 and references therein). From the theo¬ 
retical side, the equilibrium nuclear shapes in this mass 
region have been shown to suffer rapid changes as a func¬ 
tion of the number of nucleons with competing spherical, 
axially symmetric prolate and oblate, and triaxial shapes 
at close energies. Both relativistic ii and nonrelativis- 
tic approaches agree in the general description of the 

nuclear structural evolution in this mass region, which is 

S ported experimentally by spectroscopic studies [13- 
2+ lifetime measurements Il3l - ll5l| and quadrupole 
moments for rotational bands [15l| . as well as by laser 
spectroscopy measurements M- 

However, the nuclear structure richness is not the only 
attractive feature characterizing these nuclei. Another 
remarkable property of nuclei in this mass region is that 
they are involved in the astrophysical rapid neutron cap¬ 
ture process (r process), which is considered to be one 
of the main nucleosynthesis mechanisms leading to the 
production of heavy neutron-rich nuclei in the universe 
[ItL [T3| . The r-process nucleosynthesis involves many 
neutron-rich unstable isotopes, whose neutron capture 
rates, masses, and /3-decay half-lives {T 1 / 2 ) are crucial 
quantities to understand the possible r-process paths, the 
isotopic abundances, and the time scales of the process 
[ISl - l^ . Although much progress has been done measur- 
ing masses (see for example the Jyvaskyla mass database 
1^ ) and half-lives [22 h^ . most of the nuclear properties 
of relevance for the r process are experimentally unknown 
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due to their extremely low production yields in the lab¬ 
oratory. Therefore, reliable nuclear physics models are 
required to simulate properly the r process. 

The quasiparticle random-phase approximation 
(QRPA) is considered a well suited model to de¬ 
scribe medium-mass open-shell nuclear properties and 
specifically /3-decay properties. QRPA calculations for 
neutron-rich nuclei have been carried out within different 
spherical formalisms, such as Hartree-Fock-Bogoliubov 
(HFB) [ 231 , continuum QRPA with density function¬ 
als [ 231 , and relativistic mean field approaches [l^l . 
However, the mass region we are dealing with requires 
nuclear deformation as a relevant degree of freedom 
to characterize the nuclear structure involved in the 
calculation of the /3-strength functions. The deformed 
QRPA formalism was developed in Refs. (23 - [^ . using 
phenomenological mean fields. A Tamm-Dancoff ap¬ 
proximation with Sk3 interaction was also implemented 
in Ref. [s^l ■ More recently, deformed QRPA calculations 
using deformed Woods-Saxon potentials and realistic 
CD-Bonn residual forces have been performed in [sl, [111 . 
First-forbidden transitions were also considered in those 
references, showing that their effect in this mass region 
can be neglected. Various self-consistent deformed 
QRPA calculations to describe the /3-decay properties, 
either with Skyrme (33| or Gogny [^ interactions are 
also available in the literature. 

In Refs. [13, nil the decay properties of neutron- 
rich Zr and Mo isotopes were studied within a deformed 
proton-neutron QRPA based on a self-consistent Hartree- 
Fock (HF) mean field formalism with Skyrme interactions 
and pairing correlations in BCS approximation. Resid¬ 
ual spin-isospin interactions were also included in the 
particle-hole and particle-particle channels [13, Ii3|- In 
this work this study is extended to the neighboring re¬ 
gions including even and odd-A neutron-rich Ge, Se, Kr, 
Sr, Ru, and Pd isotopes. These calculations are timely 
because they address a mass region which is at the bor- 
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derline of present experimental capabilities for measur¬ 
ing half-lives at MSU and RIKEN [22l - [^ . In addition, 
theoretical calculations can be tested with the available 
experimental information on half-lives providing simul¬ 
taneously predictions for the underlying Gamow-Teller 
strength distributions and for the half-lives of more exotic 
nuclei not yet measured. Finally, this more comprehen¬ 
sive study allows one to judge better the extent to which 
the method is able to describe the decay properties of nu¬ 
clei in a wider mass region that includes spherical, well 
deformed, and weakly deformed transitional isotopes, as 
well as isotopes exhibiting shape coexistence. Therefore, 
the theoretical method will be tested over a rich set of dif¬ 
ferent nuclear structures that will reveal the limitations 
of the model. 

The paper is organized as follows. In Sec. HIl a review 
of the theoretical formalism used is introduced. Sec¬ 
tion imi contains the results obtained for the potential 
energy curves (PEC), Gamow-Teller (GT) strength dis¬ 
tributions, and /3-decay half-lives, which are compared 
with the experimental data. Section IV summarizes the 
main conclusions. 


II. THEORETICAL EORMALISM 

A summary of the theoretical framework used in this 
paper to describe the /3-decay properties in neutron-rich 
isotopes is shown in this section. More details of the 
formalism can be found elsewhere [s^ . The method 

starts from a self-consistent calculation of the mean field 
by means of a deformed Skyrme Hartree-Fock procedure 
with pairing correlations in BCS approximation. Single¬ 
particle energies, wave functions, and occupation ampli¬ 
tudes are generated from this mean field. The Skyrme 
interaction SLy4 [ilj is used as a representative of mod¬ 
ern Skyrme forces. It has been very successful describing 
nuclear properties all along the nuclear chart and has 
been extensively studied p 0, II 3 ■ 

The solution of the HE equation, assuming time rever¬ 
sal and axial symmetry, is found by using the formalism 
developed in Ref. . The single-particle wave functions 
are expanded in terms of the eigenstates of an axially 
symmetric harmonic oscillator in cylindrical coordinates, 
using twelve major shells. The pairing gap parameters 
for protons and neutrons in the BCS approximation are 
determined phenomenologically from the odd-even mass 
differences In a further step, constrained HE cal¬ 
culations with a quadratic constraint are performed to 
construct the PECs, analyzing the nuclear binding ener¬ 
gies in terms of the quadrupole deformation parameter 
/3. Calculations for GT strengths are performed subse¬ 
quently for the various minima in the energy curves in¬ 
dicating the equilibrium shapes of each nucleus. Since 
decays connecting different shapes are disfavored, simi¬ 
lar shapes are assumed for the ground state of the parent 
nucleus and for all populated states in the daughter nu¬ 
cleus. The validity of this assumption was discussed for 


example in Refs. [30l|. 

To describe GT transitions, a separable spin-isospin 
residual interaction in the particle-hole (ph) and particle- 
particle (pp) channels is added to the mean field and 
treated in a deformed proton-neutron QRPA [28l - [3^ 1^ . 
Idol . IdSjl . An optimum set of coupling strengths could be 
chosen following a case by case fitting procedure and one 
will finally get different answers depending on the nu¬ 
cleus, shape, and Skyrme force. However, since the pur¬ 
pose here is to test the ability of QRPA to account for 
the GT strength distributions in this mass region with 
as few free parameters as possible, the same coupling 
strengths are used for all the nuclei considered in this 
paper, which are taken from previous works (3^.13^. We 
use x^T — 0-1® MeV and k^j. = 0.03 MeV for the resid¬ 
ual interaction in the ph and pp channels, respectively. 

The sensitivity of the GT strength distributions to 
the various ingredients contributing to the deformed 
QRPA calculations, namely to the nucleon-nucleon effec¬ 
tive force, to deformation, to pairing correlations, and to 
residual interactions, have been investigated in the past 
[ 3 ^ HO, IdOl - li^ . In this work the most reasonable choices 
found in those references are used. Summarizing the var¬ 
ious sensitivities, the conclusion is that the main features 
of the GT strength distributions are in general very ro¬ 
bust against the Skyrme force used, showing some more 
sensitivity in the spherical cases, where the location of the 
single-particle energies is more critical to determine the 
excitation energies of the GT transitions. Deformation 
has been shown to be an important issue to describe the 
profiles of the GT strength distributions. First, because 
the degeneracy of the spherical shells is broken making 
the GT strength distributions more fragmented than the 
corresponding spherical ones. Secondly, because the en¬ 
ergy levels of deformed orbitals cross each other in a way 
that depends on the magnitude of the quadrupole defor¬ 
mation as well as on the oblate or prolate character. This 
level crossing may lead in some instances to sizable differ¬ 
ences in the GT profiles, a fact that has been exploited to 
learn about the nuclear shape from the measured /3-decay 
pattern (49l - l5i| . Pairing correlations are also important 
to describe nuclei out of closed shells. Their influence on 
the GT profiles was studied in Ref. [i^ , concluding that 
the main effect is to decrease slightly the strength at low 
energies and to create new transitions, mainly at high 
energies, that are forbidden in the absence of such corre¬ 
lations. The effect of the ph and pp residual interactions 
is also well known. The repulsive ph interaction redis¬ 
tributes the GT strength by shifting it to higher excita¬ 
tion energies causing a displacement of the GT resonance. 
It also reduces somewhat the total strength. The attrac¬ 
tive pp interaction moves the strength to lower energies. 
Its effect on the GT resonance is in general negligible, 
but nevertheless, the changes induced in the low-energy 
region are of great relevance in the calculation of the /3- 
decay half-lives, which are only sensitive to the strength 
contained in the energy region below the Q-energy win¬ 
dow. 
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The GT transition amplitudes in the intrinsic frame 
connecting the ground state |0“'') of an even-even nucleus 
to one phonon states in the daughter nucleus \uik) (K = 
0 ,1) are found to be 


( 1 ) 

where 

^ ) , ( 2 ) 

MX- = Y, {q^^XXX + , ( 3 ) 

with 

qizu — , qTTV — , ( 4 ) 

in terms of the occupation amplitudes for neutrons and 
protons ,r = 1 ~ the matrix elements of 

the spin operator, = {v\aK\T^), connecting proton 
and neutron single-particle states, as they come out from 
the HF-I-BCS calculation. X^- and Y^’- are the forward 
and backward amplitudes of the QRPA phonon operator, 
respectively. 

Once the intrinsic amplitudes in Eq. © are calculated, 
the GT strength B^{GT^) in the laboratory system for 
a transition /^^^(O+O) —>■ can be obtained 

as 


B^(Gr±) 


Y [(w/f=o \crot^ \ 0)^ S{UJK=0 - w) 

UJK 


+2{uJK=i\crit^\0y S{ujK=i - uj) ,( 5 ) 


in [5^/47r] units. To obtain this expression, the initial 
and final states in the laboratory frame have been ex¬ 
pressed in terms of the intrinsic states using the Bohr- 
Mottelson factorization [b^ . 

The s pec ific treatment of odd-B systems has been de¬ 
scribed (sTI . by considering two types of GT contri¬ 
butions. One type is due to phonon excitations in which 
the odd nucleon acts only as a spectator. The transition 
amplitudes in the intrinsic frame are in this case basically 
the same as in the even-even case, but with the blocked 
spectator excluded from the calculation. The other type 
of transitions involves the odd nucleon and is treated per- 
turbatively by taking into account phonon correlations to 
first order in the quasiparticle transitions. The excitation 
energies of the GT states with respect to the ground state 
in the daughter nuclei have been discussed in Ref. [i^ for 
both types of transitions in terms of the QRPA phonon 
energy and the quasiparticle energies. 
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FIG. 1: Potential energy curves for even-even neutron-rich 
Ge, Se, Kr, Sr, Zr, Mo, Ru, and Pd isotopes obtained from 
constrained HF-I-BCS calculations with the Skyrme force 
SLy4. 










































































4 



0.5. 
0.4 - (a) 
0.3 


—'-1-'-1-'-1-'-1~ 

3 / 2 + 5 / 2 - 5 / 2 + 

0 . 2 ^ ^ 
CQ. 0.1-®^ f rKr)=5/2(+) f fKr)=l/2(+) J 

exp exp 

.01 I C/”Kr)=l/2+ f^^rKr)=(3/2+) 

-0.2 — 3^+^'®-^. 


7/2+ 




11/2- 3/2+ 


1/2+ 5/2+ 9/2- — 

I ' I ' I ' I 



prolate 
♦ oblate 

o g-s- 

• exp 


A 


96 98 100 102 104 

A 


FIG. 2: (Color online) Isotopic evolution of the quadrupole 
deformation parameter P (a) and charge radius (b) corre¬ 
sponding to the energy minima obtained from the Skyrme 
interaction SLy4 for Ge isotopes. Ground state results are 
encircled. 


FIG. 4: (Color online) Same as in Fi g. [H but for Kr isotopes. 
Experimental charge radii are from [^ . 
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FIG. 3: (Color online) Same as in Fig. [2] but for Se isotopes. 
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FIG. 5: (Color online) Same as in F ig. [2l but for Sr isotopes. 
Experimental charge radii are from [^ . 
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FIG. 6: (Color online) Same as in Fig. [2] but for Ru isotopes. 
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FIG. 8: (Color online) QRPA-SLy4 Gamow-Teller strength 
distributions for Ge isotopes as a function of the excitation 
energy in the daughter nucleus. The calculations correspond 
to the various equilibrium deformations found in the PECs. 
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FIG. 7: (Color online) Same as in Fig. [2j but for Pd isotopes. 
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FIG. 9: (Color online) Same as in Fig. [8] but for Se isotopes. 
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FIG. 10: (Color online) Same as in Fig. |8] but for Kr isotopes. 
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FIG. 12: (Color online) Same as in Fig. [SI but for Ru isotopes. 
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FIG. 11: (Color online) Same as in Fig. E but for Sr isotopes. 13: (Color online) Same as in Fig. H but for Pd isotopes. 
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The /?-decay half-life is obtained by summing all the 
allowed transition strengths to states in the daughter nu¬ 
cleus with excitation energies lying below the correspond¬ 
ing Q-energy, Qp = Qp- = M{A, Z)-M{A, Z+ l)-me, 
written in terms of the nuclear masses M{A, Z) and the 
electron mass (me), and weighted with the phase space 
factors f{Z,Qp — E ex ); 

E f(Z,Qp-Ee,)B(GT,Ee,), 

0<Eex<Q^ 

( 6 ) 

with D = 6200 s and (gA/9v)ee = Q.77(gA/gv)tree, 
where 0.77 is a standard quenching factor. The same 
quenching factor is included in all the figures shown 
later for the GT strength distributions. The bare re¬ 
sults can be recovered by scaling the results in this paper 
for B(GT) and T 1/2 with the square of this quenching 
factor. 

The Fermi integral f(Z, Qp — Eex) is computed numer¬ 
ically for each value of the energy inclu ding screening and 
finite size effects, as explained in Ref. |53j |. 

rWo 

f (Z,Wo)= J pW{Wo-W)^X^(Z,W)dW, (7) 

with 

A±(Z, IT) = 2(1 + , 

^ ( 8 ) 

where 7 = \/l — (aZ)"^ ; y = aZWjp ; a is the fine struc¬ 
ture constant and R the nuclear radius. W is the total 
energy of the /3 particle, Wp is the total energy available 
in TOeC^ units, and p = VlT^ — 1 is the momentum in 
TOeC units. 


III. RESULTS AND DISCUSSION 

In this section I present first the PECs in the isotopic 
chains studied. Quadrupole deformation parameters as 
well as charge r.m.s. radii (re) are analyzed as a function 
of the mass number. Then, energy distributions of the 
GT strength corresponding to the local minima in the 
PECs are calculated. Finally, half-lives are evaluated and 
compared with the experiment. 

A. Structural isotopic evolution 

In Fig. [T]the PECs, i.e., the energies relative to that 
of the ground state, are plotted as a function of the 
quadrupole deformation /3 for the neutron-rich Ge, Se, 
Kr, Sr, Zr, Mo, Ru, and Pd isotopes. The results corre¬ 
spond to the SLy4 interaction. The isotopes covered in 


this study include middle-shell nuclei with proton num¬ 
bers between shell closures Z = 28 and Z = 50, namely 
Z = 32, 34,36, 38,40,42,44,46 and neutron numbers be¬ 
tween shell closures N = 50 (as in ®^Ge) and N = 82 (as 
in the heaviest ^^®Pd). 

In most of the isotopic chains one can see the appear¬ 
ance of several equilibrium nuclear shapes, whose relative 
energies change with the number of neutrons. In Ge iso¬ 
topes, prolate shapes that are ground states in the lighter 
isotopes are found with the only exception of ®^Ge, where 
a spherical shape is found in accordance with its N = 50 
semi-magic character. At N = 58,60 (®°’®^Ge) oblate 
and prolate shapes are practically degenerate in energy 
and oblate shapes become ground states for heavier iso¬ 
topes. The case of Se isotopes is similar with oblate and 
prolate minima all along the isotopic chain. The lighter 
(heavier) isotopes have prolate (oblate) ground states 
with transitional isotopes around N = 58, 60 (®^’®^Se). In 
this case the energy barriers are more pronounced than 
in the case of Ge isotopes. Kr isotopes show compet¬ 
ing shapes in the lighter isotopes that become oblate at 
N = 58,60 (®^’®®Kr) and then turn into prolate shapes 
beyond ®®Kr. In the heavier isotopes, as in the case of 
Se isotopes, shape coexistence is found with very well de¬ 
veloped oblate and prolate minima separated with high 
energy barriers. Sr isotopes show a transition from oblate 
at = 58 (®®Sr) to prolate at A = 60 (®®Sr) with a two 
minima structure for heavier isotopes. The cases of Zr 
and Mo isotopes were discussed in Refs. [13, [1^. Both 
oblate and prolate minima are observed in the lighter iso¬ 
topes of Zr and Mo with prolate ground states. Whereas 
the prolate shape remains ground state in most of the 
heavier Zr isotopes, oblate shapes are lower in energy for 
the heavier Mo isotopes. Finally, Ru and Pd isotopes 
show oblate and prolate minima in the lighter isotopes 
and a gradual transition into spherical shapes as one ap¬ 
proaches the shell closure at N = 82. 

In summary, a large diversity of nuclear structures are 
found in this mass region, from spherical to well deformed 
shapes, passing through soft transitional nuclei and even 
possible shape-coexistence structures. This rich variety 
of shapes represents a challenge to any theoretical model 
trying to describe them in a unified manner. In the next 
subsections the results are compared with the available 
experimental data, which are restricted at present to ,5- 
decay half-lives. Then, the theoretical approach will be 
tested against this information and the limitations of the 
model will be established. 

These results are in qualitative agreement with sim¬ 
ilar calculations obtained in this mass region from dif¬ 
ferent theoretical approaches, including macroscopic- 
microscopic methods based on liquid drop models with 
shell corrections [13, [ 111 , relativistic mean fields [HI] , as 
well as nonrelativistic calculations with Skyrme [ 1 ] and 
Gogny interactions. Thus, a consistent theoret¬ 

ical description emerges, which is supported W the still 
scarce experimental information available [1, [lol - [l6l IHsjl . 

The isotopic evolution can be better appreciated in 






Figs. [na where quadrupole deformations /3 (a) and 
r.m.s. charge radii (b) of the various energy minima 
are plotted as a function of the mass number A. The 
deformation corresponding to the ground state for each 
isotope is encircled. Also shown in these figures for odd- 
A isotopes, are the spin and parity (J^) of the different 
shapes and the experimental assignments [d^. The ex¬ 
perimental assignments based on systematics estimated 
from trends in neighboring nuclides have not been in¬ 
cluded. 

In Fig. [^for Ge isotopes one can see clearly the shape 
transition at A = 90 — 92 [N = 58 — 60) from prolate 
shapes with /3 Ri 0.2 to oblate shapes with /3 r; —0 .2. 
Charge r.m.s. radii have not been measured in these 
isotopes, but it is expected from these calculations a very 
smooth behavior given that the magnitude of /3 in the 
prolate and oblate sectors are very similar. Fig. [3] shows 
the analogous results for Se isotopes. In this case one 
can see the transition from prolate (/3 r; 0.2) to oblate 
{j3 K. —0.2) at A = 92 (A^ = 58). The prolate shape 
grows in the heavier isotopes (/3 Ri 0.3), but they are 
never ground states and then, the expected increasing in 
the charge radii is smooth. Kr isotopes in Fig. |4] show 
first a shape transition from prolate ~ 0.15) to oblate 
(/3 r; —0.25) and a subsequent transition from oblate 
to prolate (/3 r: 0.35) shapes. The radii are sensitive to 
this transitions, although the measured radii (s^ seem to 
favored prolate shapes in the lighter isotopes. Sr isotopes 
in Fig. [5] show a clear transition from oblate to strong 
prolate (/3 r; 0.4) deformations at A = 96 — 98 {N = 
58 — 60). This shape transition is well correlated with 
the trend change observed in the charge radii that shows 
a sizable jump between ®®Sr and ®®Sr both theoretical 
and experimentally . In the case of Ru (Pd) isotopes 
shown in Fig. IHGl), one can see a smooth transition from 
deformed oblate (prolate) solutions in the lighter isotopes 
to spherical shapes in the heavier ones. This change is 
felt in the trends of the radii, but no experiments are yet 
available to compare with. 

Spins and parities in odd-A isotopes can be compared 
with their experimental assignments. In the Ge isotopes 
the calculations agree reasonably well with the assign¬ 
ments taking into account that oblate and prolate shapes 
are very close in energy and that a I/2''' isomer is ob¬ 
served experimentally in ®^Ge at 248 keV. In the lighter 
Se isotopes, 1/2+ and 3/2+ states are obtained, whereas 
experimental assignments are (5/2+). In both isotopes, 
5/2+ states very close in energy to the ground states 
are also obtained, although somewhat above. Similarly, 
in the lighter Kr isotopes the experimental assignments 
are obtained very close in energy to the ground states, 
although slightly above. On the other hand, a 7/2+ iso¬ 
mer is experimentally observed in ®^Kr at 355 keV that 
corresponds to the ground state here. Sr isotopes ex¬ 
hibit a nice agreement. The measured spin and parities 
of ground states in 95.99.ioigj. correspond to the prolate 
calculations. A (7/2+) state is also observed experimen¬ 
tally in ®®Sr at 56 keV. In ®^Sr the observed 1/2+ ground 


state appears as an excited state. It is also worth noting 
that the prolate ground state (3/2“) for this isotope is 
observed experimentally at 645 keV. In the case of Ru iso¬ 
topes the measured are difficult to reproduce. They 
are found in the calculations, but not as ground states. 
On the other hand, the negative parity 7/2“ states found 
in the calculations are also seen experimentally at low en¬ 
ergies. In particular, an isomeric state (7/2“) at an un¬ 
determined energy has been seen in ^^^Ru. Finally, in Pd 
isotopes the negative parity isomers, which are oblate in 
this description, are reproduced in the calculations, but 
not the ground states. 
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FIG. 14: (Color online) QRPA-SLy4 accumnlated GT 
strengths in Ge isotopes calcnlated for the various eqnilib- 
rium shapes. Qp and Sn energies are shown by solid and 
dashed vertical arrows, respectively. 


B. Gamow-Teller strength distributions 

In the next figures, the energy distributions of the GT 
strength corresponding to the various deformed equilib¬ 
rium shapes are shown for each isotopic chain. The re¬ 
sults are obtained from QRPA with the force SLy4 with 
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FIG. 15: (Color online) Same as in Fig. 1141 but for Se iso¬ 
topes. 


FIG. 16: (Color online) Same as in Fig. [Til but for Kr iso¬ 
topes. 


pairing correlations and with residual interactions with 
the parameters written in Sec. im The GT strength is 
plotted versus the excitation energy of the daughter nu¬ 
cleus with a quenching factor 0.77. Zr and Mo isotopes 
were already studied in Refs. p, m and are not re¬ 
peated here. 

Figures contain the results for Ge, Se, Kr, Sr, 
Ru, and Pd isotopes. The energy distributions of the in¬ 
dividual GT strengths corresponding to the ground state 
shapes are shown, together with continuous distributions 
for the ground state shapes as well as for the other possi¬ 
ble shapes, obtained by folding the strength with 1 MeV 
width Breit-Wigner functions. Qp values are shown with 
vertical arrows. In both cases, even and odd isotopes, 
the Qp values increase with the number of neutrons in 
each isotopic chain and the values in the odd-A isotopes 
{Z, N + 1) are about 2-3 MeV larger than the values 
in the neighbor even-even isotopes {Z,N). The general 
structure of the GT distributions is characterized by the 
existence of a GT resonance, which is placed at increasing 
excitation energy as the number of neutrons N increases 
in a given isotopic chain. The total GT strength also 


increases with TV, as it is expected to fulfill the Ikeda 
sum rule. The various shapes produce quite similar GT 
strength distributions on a global scale. Nevertheless, the 
small differences among the various shapes at the low en¬ 
ergy tails (below the Qp) oi the GT strength distributions 
that can be appreciated because of the logarithmic scale, 
lead to sizable effects in the /3-decay half-lives. 

Unfortunately, comparison with experiment is still not 
possible for the GT strength distributions, the measured 
half-lives will be compared to the calculations in the next 
subsection. Gomparison with calculated GT distribu¬ 
tions from other theoretical approaches is also restricted 
to the few cases where these results have been published 
In Refs. [11,163 the authors performed QRPA 
calculations with deformed Woods-Saxon potentials and 
realistic GD-Bonn residual forces using the G-matrix for¬ 
malism and compared these results with the results ob¬ 
tained from separable forces. While in Ref. [Ill the 
comparison between the results obtained from realistic 
or separable residual interactions is restricted to the half- 
lives, in Ref. [60l| the authors compared those results in 
the context of two-neutrino double-beta decay, conclud- 
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FIG. 17: (Color online) Same as in Fig. 1141 but for Sr iso¬ 
topes. 


FIG. 18: (Color online) Same as in Fig. 1141 but for Ru iso¬ 
topes. 


ing that both approaches, realistic and separably lead 
to similar results. On the other hand, in Ref. the 
Skyrme force SLy4 was used to generate the mean field 
as it is done in this work. The residual interaction in 
the ph channel was self-consistently introduced and not 
reduced to a separable form. Finally the pp residual in¬ 
teraction was written as a contact force with a coupling 
strength fitted to reproduce the half-life in ^™Zr. The 
GT strength distributions in neutron-rich Zr isotopes ob¬ 
tained from this approach were compared with the cor¬ 
responding distributions obtained with separable forces 
in Figs. 5-6 in Ref. [s^. From this comparison one can 
conclude that in many aspects the main characteristics of 
the consistent force are maintained by a separable force 
with a much lower computational cost. The compari¬ 
son of the half-lives shows also a remarkable agreement 
between both approaches. 

In the next figures, Figs. [TiHTOl one can see in more 
detail the accumulated GT strength in the energy re¬ 
gion below the corresponding Q/S energy of each isotope, 
which is the relevant energy range for the calculation of 
the half-lives. The vertical solid (dashed) arrows show 


the Qp (Sn) energies, taken from experiment (43 |. In 
these figures the sensitivity of these distributions to de¬ 
formation can be appreciated and one can understand 
that measurements of the GT strength distribution from 
/3-decay can be, in particular cases, an additional source 
of information about the nuclear deformation, as it was 
shown in Refs, The GT strength distribution in 

odd-A isotopes is found to be displaced to higher energies 
(typically about 2-3 MeV) with respect to the even-even 
case. The shift corresponds roughly to the breaking of a 
neutron pair and therefore it amounts to about twice the 
neutron pairing gap. Below this energy only transitions 
involving the odd nucleon are possible. 

The energy distribution of the GT strength is funda¬ 
mental to constrain the underlying nuclear structure. For 
a theoretical model, it represents a more demanding test 
than just reproducing half-lives or total GT strengths 
that are integral quantities obtained from these strength 
distributions properly weighted with phase factors (see 
Eq. ini). These quantities might be reproduced even with 
wrong strength distributions. This is of especial impor¬ 
tance in astrophysical scenarios of high densities and tern- 
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FIG. 19: (Color online) Same as in Fig. [141 but for Pd iso¬ 
topes. 


peratures that cannot be reproduced in the laboratory. 
Given that the phase factors in the stellar medium are dif¬ 
ferent from those in the laboratory, the stellar half-lives 
become dependent on the electron distribution in the 
stellar plasma that eventually may block the /3-particle 
emission [^. Therefore, to describe properly the decay 
rates under extreme conditions of density and tempera¬ 
ture, it is not sufficient to reproduce the half-lives in the 
laboratory. One needs, in addition, to have a reliable 
description of the GT strength distributions [bl, . 


C. Beta-decay Half-lives 

The calculation of the half-lives in Eq. (0) involves 
knowledge of the GT strength distribution and of the /3 
energies (Qp — Ee^), which are evaluated by using Qp 
values obtained from the mass differences between par¬ 
ent and daughter nuclei obtained from SLy4 with a zero- 
range pairing force and Lipkin-Nogami obtained from the 
code HFBTHO [H. 

In Figs. HOHSnithe measured /3-decay half-lives (solid 
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FIG. 20: (Color online) Measured /3-decay half-lives for Ge 
isotopes compared to theoretical QRPA-SLy4 results calcu¬ 
lated from different shapes. Circles are experimental values 
(open circles are experimental values from systematics) [43 |. 


dots, open dots stand for experimental values from sys¬ 
tematics) [1^ are compared with the theoretical re¬ 
sults obtained with the prolate, oblate, and spherical 
equilibrium shapes, for the various isotopic chains. In 
Fig. [20] one can see the half-lives for Ge isotopes. The 
lighter isotopes are not well reproduced, being largely 
overestimated. This point will be discussed later. The 
half-lives obtained from oblate shapes are larger than the 
corresponding prolate ones. This feature is correlated 
with the GT strength contained below the Qp energy in 
Figs. 151 and im Prolate shapes, which are closer to exper¬ 
iment, are also the ground states in this range of masses 
according to the calculations (see Fig. [5]). For heavier 
isotopes, the half-lives for oblate and prolate shapes are 
very similar. In the case of Se isotopes in Fig. |2I1 the 
calculations also overestimate the half-lives of the lighter 
isotopes, but the agreement with experiment is in this 
case much better. In the middle region the experimental 
half-lives, which are taken from systematics, are reason¬ 
ably well reproduced. The half-lives of heavier isotopes 
exhibit a rather flat behavior. Half-lives of Kr isotopes 
are shown in Fig. [22l As in the previous figures, the 
half-lives from the oblate shapes are larger than the pro¬ 
late ones in the lighter Kr isotopes, but the situation is 
reversed at ®"'^Kr. This is again nicely correlated with the 
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GT strength at low excitation energies shown in Fig. [161 
In general, the half-lives in the middle region are well de¬ 
scribed. This is also true for Sr and Ru isotopes in Figs. 
|23]and|24l respectively, where the trends observed exper¬ 
imentally are well reproduced, except for the lighter Sr 
isotopes that are clearly underestimated and the heavier 
Ru isotopes, where the data from systematics fall down 
faster than the calculations. Finally, in the case of Pd 
isotopes, shown in Fig. [551 the calculations underesti¬ 
mate (overestimate) the measured half-lives in the lighter 
(heavier) isotopes. 

All in all, the agreement with experiment is reason¬ 
able, especially in the middle regions. These regions con¬ 
tain in general well deformed nuclei, where the present 
approach is more suitable. On the other hand, weakly 
deformed transitional isotopes, such as light Ge and Se 
isotopes and heavy Ru and Pd isotopes are not so well 
described. Furthermore, in the light isotopes of all the 
isotopic chains, which are closer to the valley of stability, 
the half-lives are larger because of the small Q /3 energies 
involved. In these cases the half-lives are determined ex¬ 
clusively by the very low energy tail of the GT strength 
distribution contained in the narrow window below Qp. 
Therefore, tiny variations in the description of the GT 
strength distribution in the low-lying energy region can 
drive sizable effects in the half-lives. Of course it is also 
important to describe the half-lives of the long-lived iso¬ 


topes, but their significance to constrain the GT strength 
distribution is minor since the half-lives are insensitive to 
most of this distribution. 

Half-lives for neutron-rich Kr, Sr, Zr, and Mo isotopes 
calculated from self-consistent deformed QRPA calcula¬ 
tions with the Gogny DIM interaction and experimental 
values of < 5/3 agree with the results in this work within 
the uncertainties of the calculations. The agreement is 
also very reasonable between the calculated half-lives and 
those obtained from deformed QRPA calculations using 
deformed Woods-Saxon potentials to generate the mean 
field and complemented with realistic CD-Bonn residual 
forces [H, I 33 . T he agreement is also good with the re¬ 
sults in Ref. using the Skyrme force SLy4 with con¬ 
sistent residual interactions in the ph channel as men¬ 
tioned earlier. Fig. 7 in that reference displays this com¬ 
parison. 

It is also worth noticing that the worst agreement with 
experiment occurs in the light Ge isotopes, as well as 
in heavy Pd isotopes. In these cases the calculations 
overestimate the experiment leaving room for contribu¬ 
tions coming from first forbidden (FF) transitions. One 
can understand from simple qualitative arguments that 
the role of FF transitions is expected to be more impor¬ 
tant in lighter Ge and in Pd isotopes. Thus, for 32 Ge, 
3480 , 36 Kr, and 38 Sr isotopes, the last occupied proton 
orbitals come basically from the 2 ^ 3 / 2 , 1/572 and 2 pi /2 
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negative-parity spherical shells. On the other hand the 
neutrons in isotopes occupy orbitals belonging to 

the 1 ( 79 / 2 , 2 ^ 5/2 and 157/2 positive-parity spherical shells. 
Therefore, in the /3-decay one neutron in a positive-parity 
state is transformed into a proton that would sit in a 
negative-parity state, thus suppressing GT and favoring 
FF transitions in the low-lying transitions. This is par¬ 
ticularly true for the lighter Ge isotopes. In the heavier 
ones, other neutron states with negative parity (I/ 111 / 2 ) 
have to be considered because of deformation effects. The 
same argument can be applied to the lighter Se, Kr, and 
Sr isotopes, but in these cases proton states from positive 
parity (I 59 / 2 ) are closer in energy and would participate 
in the decay favoring GT transitions. The situation is 
different in the case of Ru and Pd isotopes. Now the 
available proton states for the decay are of positive parity 
( 159 / 2 ), whereas most of the last occupied neutrons be¬ 
long to negative-parity states (I/ 111 / 2 ), thus favoring FF 
transitions. According to calculations [^ . [s^ of the FF 
transitions in this mass region, minor effects are expected 
from them. Nevertheless, it would be very interesting in 
the future to study systematically the FF contributions 
in all the isotopes in this mass region. 

Another feature observed in the present calculations is 
the existence of some odd-even staggering effect in the 
calculated half-lives, which is not observed experimen¬ 
tally. This effect is particularly evident in Ru and Pd 


isotopes. There are not many calculations involving si¬ 
multaneously even-even and odd-A isotopes, but some of 
them exhibit some sort of staggering effect as well [sj. 
The appearance of this effect in the half-lives suggests 
some deficiency in the model that might be related to 
the determination of ground-state energies in the odd-A 
systems [h^. Unfortunately, there are more sources of 
uncertainty related to the odd-A systems that should be 
considered as well [g^, such as the spin and parity as¬ 
signments, the blocking procedures or the treatment of 
the Iqp excitations involving the odd nucleon. This is¬ 
sue will be the subject of a future investigation in this 
direction. 

It is also interesting to look for the simultaneous ap¬ 
pearance of structural effects that eventually can appear 
in different observables. One example can be seen in the 
evolution of the experimental half-lives with the number 
of neutrons in the isotopic chains. At some points one 
observes discontinuities in the general trends of behavior, 
such as in the mass regions ®°’®^Se, ®^’®‘*Kr, ®®’®®Sr, and 
ii®’120Ru. These experimental findings on the half-lives 
are correlated with the shape transitions in Figs. |3H6]pre- 
dicted in the model. One cannot state firmly that these 
sharp changes in the behavior of the half-lives are signa¬ 
tures of shape transitions, but certainly this correlation 
cannot be discarded given that a change of the deforma¬ 
tion in the nuclear system involves a structural change 
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to whom the half-lives are also sensitive. 

Finally, the impact of deformation on the decay prop¬ 
erties can be better appreciated in a systematic compar¬ 
ison of the half-lives calculated with both the spherical 
approximation and the deformation that corresponds to 
the minimum of the PEC for each isotope. Then, Fig. 
[26l shows the ratios of the calculated and experimental 
half-lives for two sets of data corresponding to a spher¬ 
ical calculation (open dots) and to a deformed calcula¬ 
tion (solid dots) at the self-consistent deformation that 
gives the minimum of the PECs. These ratios are plotted 
as a function of the experimental half-lives (a) and as a 
function of the quadrupole deformation at the minimum 
of the PECs (b). To increase the size of the sample, 
besides the isotopes considered in this work with mea¬ 
sured half-lives, I have also included the set of Zr and 
Mo neutron-rich isotopes studied in Ref. [s^ with mea¬ 
sured half-lives. In the upper panel of Eig. [26] (a) one 
can see how deformation improves the description of the 
half-lives. Practically all the full black dots are contained 
within the horizontal lines dehning the region of one or¬ 
der of magnitude agreement. On the other hand, the 
results from the spherical calculation are more spread 
out with larger discrepancy with experiment. One can 
also see that the results are better in both spherical and 
deformed calculations for shorter half-lives, whereas the 
results for larger half-lives show sizable deviations. The 


FIG. 26: (Color online) Ratio of calculated to experimental 
d-decay half-lives for two sets of calculations, with the spheri¬ 
cal approximation (open dots) and with the deformation that 
corresponds to the minimum of the PECs (solid dots). The 
ratios are plotted as a function of the experimental half-lives 
(a) and as a function of the quadrupole deformation at the 
minimum of the PECs (b). 


latter correspond to isotopes close to the valley of sta¬ 
bility with small Q^-values, where the half-lives are only 
sensitive to the small portion of the GT strength distri¬ 
bution at low excitation energies below Qp. In the lower 
panel (b) one can see the results from a different point of 
view and it can be studied whether deformation improves 
the results evenly in the whole range of deformations or 
whether its effect is stronger at large deformations. Three 
regions of accumulation of results can be distinguished. 
Two of them correspond to well deformed nuclei located 
at /3 Ri —0.2 and /3 r; 0.35. In these regions the de¬ 
formed calculations clearly improve the results from the 
spherical ones that show a tendency to underestimate the 
experiment. The other region corresponds to 0 < /3 < 0.2 
values, where nuclei are softer or transitional and the de¬ 
formed formalism should be improved. In this region the 
results are more scattered than in the well deformed re¬ 
gions, but the deformed calculations show deviations that 
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rarely exceed one order of magnitude, still representing 
an improvement over the spherical results. 

In order to have a quantitative estimation of the 
quality of the various calculations, following the anal¬ 
ysis made in Ref. [2^, the logarithms of the ratios of 
the calculated and experimental half-lives are introduced 
through the quantities 


r = logio 


Ji/2(calc) 
.^ 1/2 (exp). 


(9) 


Then, the average position of the points, Mr, the stan¬ 
dard deviation, dr, and the total error, Sr, are defined 
as 


Mr 


1 

n 



(7 'f' 


Sr 


1 

n 





1 

n 



( 10 ) 


and their corresponding factors M^^ = 10 ^'', (t):® = 10 '^'’, 
and S);° = 10^'’. The analysis of the results shown in Fig. 
[26l involving n = 81 nuclei leads to the values M,)® = 
1.105, (Tr° = 10.21, and S);*^ = 10.24 in the spherical case 
and = 0.937, = 3.09, and S^o = 3.09 in the 

deformed one, showing clearly the improvement achieved 
with the deformed formalism. 


IV. CONCLUSIONS 

A microscopic approach based on a deformed QRPA 
calculation on top of a self-consistent mean field ob¬ 
tained with the SLy4 Skyrme interaction has been used 
to study the nuclear structure and the decay properties 
of even and odd neutron-rich isotopes in the mass region 
A Ri 80 — 130. The nuclear model and interaction have 
been successfully tested in the past providing good agree¬ 
ment with the available experimental information on bulk 
properties all along the nuclear chart. Decay properties 
in different mass regions have been well reproduced as 
well. The structural isotopic evolution has been studied 
from their PECs. Depending on the isotopic chain, a 
large variety of nuclear shapes is found, including spher¬ 
ical shapes, well developed deformed shapes, and transi¬ 
tional soft shapes. Charge radii have been also investi¬ 
gated, showing the connection between a discontinuous 


behavior in the isotopic trend with a shape transition and 
comparing the results with the available measurements 
from laser spectroscopy. Then, Gamow-Teller strength 
distributions and ,5-decay half-lives have been computed 
for the equilibrium shapes. 

The isotopic evolution of the GT strength distributions 
exhibits some typical features, such as GT resonances in¬ 
creasing in energy and strength as the number of neu¬ 
trons increases. Effects of deformation are hard to see 
on a global scale, but they become apparent in the low 
excitation energy below energies, a region that deter¬ 
mines the half-lives. Half-lives have been calculated using 
Qfj energies calculated with the force SLy4. In general, 
a reasonable agreement with experiment is obtained, es¬ 
pecially in the short-lived nuclei of Ge, Se, Kr, Sr, and 
Ru isotopes. The results are comparable to other calcula¬ 
tions using different approaches for the mean field and/or 
residual interactions. Special difficulties are found to de¬ 
scribe properly the half-lives of the lighter Ge isotopes 
and the Pd isotopes. These are examples of transitional 
nuclei where the nuclear structure is more involved and 
the concept of a well defined shape might not be mean¬ 
ingful. 

A systematic comparison of the ratios of the calcu¬ 
lated and experimental half-lives has been done using 
both spherical and deformed calculations, showing that 
the inclusion of deformation improves significantly the 
description of the decay properties. 

Experimental information on the energy distribution of 
the GT strength is a valuable piece of knowledge about 
nuclear structure in this mass region. The study of these 
distributions is within the current experimental capabil¬ 
ities in the case of the lighter isotopes considered in this 
work. Here, I have presented theoretical predictions for 
them based on microscopic calculations. Similarly, mea¬ 
suring the half-lives of the heavier isotopes will be highly 
beneficial to model the r process and to constrain theoret¬ 
ical nuclear models. This possibility is also open within 
present capabilities at RIKEN. 
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